Plasmon attenuation and optical conductivity of a two-dimensional electron gas 
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In a ballistic two-dimensional electron gas, the Landau damping does not lead to plasmon atten- 
uation in a broad interval of wave vectors q ^ kp. Similarly, it does not contribute to the optical 
conductivity a{Lu,q) in a wide domain of its arguments, Ep > u> > qvp, where Ef, kp and vp 
are, respectively, the Fermi energy, wavevector and velocity of the electrons. We identify processes 
that result in the plasmon attenuation in the absence of Landau damping. These processes are: the 
excitation of two electron-hole pairs, phonon-assisted excitation of one pair, and a direct plasmon- 
phonon conversion. We evaluate the corresponding contributions to the plasmon linewidth and to 
the optical conductivity. 

PACS numbers: 73.23.-b, 72.30.-|-q, 71.45.Gm, 63.22.+m 




I. INTRODUCTION 

Because of the long-range nature of the Coulomb inter- 
action, the energy of a plasmon propagating in an elec- 
tron gas exceeds greatly the energy of an electron-hole 
excitation at the same wave vector g, provided that q is 
smaller than the Fermi wave vector kp- This remains 
true in any dimensionality, including the case of a two- 
dimensional electron gas (2DEG), for which the plasmon 
spectrum is gapless. 



(1) 



here is the Fermi velocity, k = 2me^/e* is the inverse 
screening radius, m is the effective electron mass, and 
e* is the dielectric constant of the host material. (Here- 
inafter, we use the units with % = 1.) High plasmon 
velocity, duiq/dq 2> wf, prevents plasmon from decaying 
into an electron-hole pair and makes the Landau damp- 
ing exponentially small at low temperatures. Therefore, 
the leading contribution to the plasmon attenuation of a 
purely electronic nature comes from the plasmon decay 
into two electron-hole pairs^. Two such pairs propagat- 
ing in opposite directions can carry large energy while 
having a negligible total momentum. 

In 2DEG formed in semiconductor heterostructures, 
there are additional channels for plasmon attenuation. 
The energy and momentum conservation conditions can 
also be satisfied for processes in which a plasmon cre- 
ates an electron-hole pair and a phonon, rather than two 
electron- hole pairs. We will call an attenuation due to 
such processes a phonon-assisted Landau damping since 
it involves a single electron- hole pair. The effectiveness 
of this process is increased by a large phase space for the 
emission of phonons in the bulk of a semiconductor. 

Yet another possibility of the plasmon attenuation in 
a heterostructure arises from a conversion of the plas- 
mon into an acoustic phonon via a virtual electron-hole 
pair. Such a process is feasible due to the absence of 



momentum conservation in the direction perpendicular 
to the 2DEG plane (z-direction). This makes it possible 
to satisfy energy and in-plane momentum conservation 
conditions for a process in which a plasmon is converted 
into a phonon with appropriate value qz of the momen- 



tum along z-axis: uiq 



1, where s is the sound 



velocity. Note that processes involving optical phonons 
can be safely neglected as they have large energy thresh- 
olds (typically 35 — 50 meV) . 

In the present paper, we analyze the mentioned above 
inelastic processes leading to plasmon attenuation, and 
do not consider the effect of impurities, assuming 2DEG 
sufficiently clean. In Section II the plasmon broadening 
due to plasmon-phonon conversion is discussed. The con- 
version is mediated by a virtual electron-hole pair with a 
rate independent of temperature. In subsequent Sections 
we discuss scattering processes leading to the creation of 
real electron-hole pairs in the final state: in Section III 
the probability for plasmon scattering into a phonon and 
an electron-hole pair is derived. In Section IV we evaluate 
the rate of plasmon decay into two electron-hole pairs. In 
Section V we compare the considered mechanisms with 
each other. 



II. PLASMON-PHONON CONVERSION 

The Hamiltonian of the 2DEG interacting with 
phonons has the form. 



E 



^ 2m ^ £*|x,, - Xi 

d^k Mx{^) 



(27r)3 



(cAk' 



i>j 



+ C.C. 
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where the second quantization representation for phonon 
variables and coordinate representation for electron vari- 
ables are used. The second term in Eq. Q stands 
for the electron-electron interaction and the last term 
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describes interaction of two-dimensional electrons with 
three-dimensional phonons: M\(k) is the coupling func- 
tion for the A-th phonon branch, CAk are the phonon an- 
nihilation operators, k|| denotes the in-plane component 
of the phonon wave vector k, and p is the density of the 
crystal. For simplicity, we assume an isotropic phonon 
spectrum, w^k = s\k. 

The electron-phonon interaction in semiconductor het- 
erostructures is mainly due to the deformation poten- 
tial (both Si- and GaAs-based structures) and piezo- 
electric coupling (GaAs structures). The deformation 
potentia&i4, 



M 



(3) 



couples electrons to the longitudinal (l) acoustic phonons 
only, S being the deformation potential constant. 

Piezoelectric interaction (present in GaAs heterostruc- 
tures) couples electrons both to longitudinal and trans- 
verse (t) phonona^tSii, 



Mlik) = 2[ehu) 



/r2 
fc6 



+ 8 II ' 



\ = U 



(4) 



Here we utilized the conventional notation for the cou- 
pling constant, hi^. In contrast to the deformation po- 
tential interaction Eq. Q , piezoelectric coupling does not 
vanish in the long- wavelength limit k — > 0. 

The plasmon dispersion is determined by the solution 
of the equation e{ui,q)/e* = 1 — Vqxii^,q) = 0. Here 
Vq = 27re^/(qe*) stands for the Fourier transform of the 
Coulomb interaction. The electron polarization operator 
in the random phase approximation (RPA) is given by 



x(w,g) = 



2(Pp 



{2n) 



IT] 



(5) 



where Up is the Fermi-Dirac distribution function. In the 
plasmon frequency range (w 3> qvp) the polarization op- 
erator can be approximated by xi^^ q) ~ Taq^v\l(2'KuP'\ 
leading to the plasmon spectrum, Eq. l(T]l. 



FIG. 1: The Dyson equation for the propagator of the scalar 
potential (shown by the bold wavy line) in the lowest or- 
der in the electron phonon-interaction. The thin wavy line 
stands for the propagator evaluated in the conventional RPA 
approximation, (7(cj,g) = Vq/[1 — Vqx(i.^, ?)], the zigzag line 
represents the phonon Green function, the loop stands for the 
electron polarization operator x('^,9)- 

To find the lowest-order in electron-phonon interaction 
correction to the dielectric function, one needs to solve 
the equation for the longitudinal electric field propagator, 
see Fig. 1. The resulting dielectric function is 



e(w,q)/e* 



1 



E 



dqz Mliq,qz) 



27r p[{u! + iriY 



(6) 



The plasmon spectrum lo = ujq — ijq is determined by a 
zero of the dielectric function ©. The imaginary part 
7g originates from the poles of the phonon propagator at 
those values of momentum q^ which satisfy the energy 
conservation condition, ojq — s\^Jq^ + 9z- After simple 
calculation, with the use of the condition |(7z| ~ ^q/ s\ 3> 
q, we obtain 



A 



(7) 



The deformation potential interaction yields the contri- 
bution to the plasmon linewidth. 
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Here ujg = kpsi is the characteristic phonon frequency, 
and 



rriKp 



2n^psf 



(9) 



is the dimensionless coupling constant for the deforma- 
tion potential interaction. 

The piezoelectric coupling, contrary to the deformation 
potential, remains finite in the long- wavelength limit. 
Therefore, when present, the former often leads to much 
stronger effects than the latter. However, the piezoelec- 
tric contribution into the width Q is strongly reduced 
due to the fact that the plasmons emit phonons almost 
perpendicular to the 2DEG plane, 9 « qs/ujq <C 1. Ac- 
cording to Eq. the piezoelectric coupling is stronger 
for the transverse phonons. However, the corresponding 
contribution to the width Eq. lO is still small compared 
to the deformation potential contribution Eq. by a 
factor of {s/vfY 1- 

Rate of the direct (via a virtual electron-hole pair) 
plasmon-phonon conversion given by Eq. ((TJ, describes 
the temperature-independent plasmon broadening. To 
analyze the temperature-dependent attenuation one has 
to account for the inelastic processes resulting in creation 
of electron-hole pairs in the final state. To do this we re- 
late the plasmon attenuation to the optical conductivity. 



III. PHONON-ASSISTED LANDAU DAMPING 

The propagation of a plasmon in 2DEG is accompanied 
by an oscillating electric field, whose scalar potential is 
a plane wave in the in-plane direction x and satisfies the 
equation = outside 2DEG, 



(x, z,i) = (0oe" 



(10) 



The energy (per unit area) of the electric field (jTUI) is 
given by 



W,^- I dz (V0)2 = i^ld^of = eV-Vol', (11) 
ovr / zvr 
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where the bar denotes the in-plane average. The kinetic 
energy of electrons is equal to the energy of electric field 
(the virial theorem). The total energy of a plasmon is 
therefore, W = 2We. The plasmon attenuation 7^ is 
related to the energy dissipation rate, 



dW/dt = -2jgW = -4-fqWe. 



(12) 



On the other hand dissipation is related to the real part 
of the longitudinal optical conductivity a'{uj,q), 



dW 

~dr 



(13) 



The right-hand side in this expression is the Joule heat- 
ing. Note, that Eq. H13|) can be applied to attenuation 
of plasma eigenmodes as well as the dissipation of the 
external electric field. In the latter case q(j)o should be 
understood as the amplitude of the electric field at the 
2DEG plane. Using Eqs. l(T^ and l(T^ we can relate 
the plasmon attenuation to the optical conductivity at 



2e 



f cr'(u;5,g). 



(14) 



Equation H14|l may be viewed as the Ward identity relat- 
ing xil^i^l) 8.nd a(uj,q), i.e., the polarization operators 
with scalar and vector vertices. 

To evaluate 7^, we find first dW/dt with the help of 
the perturbation theory in the interaction of electrons 
with phonons and with electric field (j)o, see Eqs. Q and 
((TU|l . We then obtain 7^ from dW/dt with the help of 
Eqs. H13|l - I|14|l . We are interested in the dissipation of 
the high-frequency electric field, lo 3> qvp- At such fre- 
quencies, the field absorption due to an excitation of a 
single electron-hole pair (Landau damping) is forbidden 
by the energy and momentum conservation conditions. 
The lowest-order absorption process, therefore, includes 
also a creation or annihilation of a phonon, see Fig. 2. 
The probability of the field absorption process is given 




p+q 





p-q 



FIG. 2: Graphic representation of the transition amplitude 
L. Wavy lines stand for the amplitude of the electric field 00, 
zigzag lines denote phonons. 

by the Fermi Golden rula^, 

dw± = 2tt\Lx\^6{^p +oJ-^p,T WAk) 



xS{p + q - p' =F ki 



d'^p'd^k 
(27r)3 ■ 



(15) 



Here the matrix element 



1 



■Cp Cp+q 



?P'-^P 



(16) 



accounts for the intermediate state in the process, and 
the signs ± correspond to the creation (annihilation) 
of a phonon, which accompanies the excitation of the 
electron-hole pair. The dielectric function in the de- 
nominator of the matrix element Eq. (jKil) accounts for 
the screening of the electron-phonon interaction by two- 
dimensional electrons. The presence of two terms in 
Eq. Hlt)|) corresponds to the two different possibilities for 
the virtual state depending on whether the electron in- 
teracts first with the plasmon or with the phonon, see 
Fig. 2. Each internal line contributes to the matrix ele- 
ment a propagator of the virtual state, {E — Ey)~^, with 
E being the total energy of particles, and Ey the energy 
of the virtual state. 

The resulting probability of the field absorption is 
given by the appropriate sum over the initial and final 
states of the electron-phonon system. 



E 

A 



2cfp_ 



np{l - np/)[(l + Nxk)dw+ + Nxkdw- 



(17) 

where Nxk is the Bose-Einstein phonon distribution, and 
factor 2 takes into account the electron spin degeneracy. 
The probability for the field emission, /g can be found 
from the detailed balance principleSi, le = /ae~"/"^. 

The energy dissipation rate is then determined by the 
two probabilities. 



dW 

'~dr 



- Q - c^(l - e-"/^)4 



(18) 



In the high-frequency domain uj ^ qvp, the relevant mo- 
mentum of an electron-hole pair p — p' and that of a 
phonon k are large compared to the plasmon momentum 
q, which allows us to simplify the matrix element Lx, 



e(l3oe*Mxi'k) (p - pQ ■ q 
e{uj,^)y/2pLUxk muj^ 



(19) 



In the last expression we also have used the momentum 
conservation condition Eq. (|15|l . Combining Eqs. I|13|l . 
(I15|l and p7|l -H19 |l . we obtain the phonon-assisted Lan- 
dau damping in the form: 



x e^£*^sinh[^] ^ 



d^k 



k^Mlik) 



(2^)3wAfe|e(^,fc||)|^ 



1 



X"(^A/c - t^,fc||) X"('^A/c fc||) 



■ sinh 



sinh [^^^^] 



sinh 



.(20) 



Here x" is the structure factor of 2DEG, i.e., the imagi- 
nary part of the polarization operator Eq. . We are in- 
terested in the optical conductivity a,t lu,T <ti Ep, where 



4 



Ep = kp/2m is the Fermi energy. It allows us to use the 
low- frequency limit of x" j 



n 



Tr\k\\\vp 



1 - f-^) 



This limiting form is valid across the entire particle-hole 
domain of excitations, except narrow regions near its 
ends, |ri|/w_F < fc|| < 2kp — \Q\/vf. 

The electron-phonon interaction is effectively screened 
at a small momentum transfer between the two subsys- 
tems, which results in a strong frequency dependence of 
the conductivity. Indeed, since only the particle-hole 
domain is relevant in the integral H2l)|l . it suffices to 
approximate the dielectric function by its static limit, 
£(0, A:||)/e* = 1-1- K^/fcj|. The characteristic momenta in 

the integral of Eq. (|2()|l are fcy ~ max{Lu,T)/ s. There- 
fore, at cl;,T ^ the screening of the electron-phonon 
interaction is strong (here we introduced the character- 
istic frequency lOk — ks, and k is the inverse screening 
radius). We consider below the plasmon attenuation in 
the limits of weak and strong screening. 

(i) In the limit of low frequency and temperature, 
uj,T <^ uj^, corresponding to the strong screening, a 
straightforward evaluation of the integral in Eq. I|2U|) 
yields for the electron-phonon interaction via deforma- 
tion potential. 



7g 



uj^, ujq> 27rr, 

A{2T:Ty, ojq < 27rr, 



(21) 



where A ~ 1.47. Similarly, for the piezoelectric coupling. 



7? 



15 • 25 cj?u;2 



ujg> 27rT, 
S(27rT)^ LUq < 27rT, 



(22) 



with the coefficient B = 0.76. Here characteristic fre- 
quencies oJk. — KSi and lUs — kpsi depend on the material 
parameters. The deformation potential coupling strength 
r]d is given in Eq. , and the dimensionless constant for 
the piezoelectric coupling is 



2'K'^psfkp 



63 

256 '^s? 256j' ^ ' 



Clearly, in the limit of low temperature and plasmon fre- 
quency, the piezoelectric part, Eq. (|22|l . prevails over the 
deformation potential contribution Eq. (|21(l . 

(ii) We turn now to the case of a weak screening. In 
typical heterostructures the inverse screening radius k is 
of the order of fc^?, and therefore a;„ ^ lOs- We concen- 
trate on relatively large frequencies or high temperatures, 
i.e., assume that uj or T exceed lu^ and lOs^ but are still 
small compared to the Fermi energy, Ep = kp/2m. 

The leading deformation potential contribution in 
Eq. (|20|l comes from phonon states with fc|| ^ kp and 
kz ^ max(a;, T)/s 3> kp^ i.e., from the phonons propa- 
gating almost normally to the 2DEG plane. In the limit 



of high temperature or frequency we find 



7, = C^K + (27rT)^], 



(24) 



Here the numerical coefficient C depends on the screening 
parameter K/2kp, 



t/2 



c = 



3 7 isin9+^r' 



(25) 



and varies from C = 7r/12 « 0.26 to C = 0.06 when 
the screening parameter increases from K./2kp = to 
n/2kp = 1. 

Unlike the interaction via deformation potential, the 
piezoelectric mechanism of the electron-phonon interac- 
tion is ineffective for phonons with |fc||| ^ k. As the 
result, the piezoelectric mechanism yields a contribution 
to 7g, which is smaller than the one given by Eq. H24|l . 
In particular, at low temperatures and high frequencies. 



< 



w, this contribution to the plasmon line width 



saturates at a value 

Iq ^vl^s, (26) 
while in the case of high temperature, T ^ LUg, we find 



2 rj. 



(27) 



independent of the parameter lo/lOs. Strictly speaking, 
the piezoelectric constants in Eqs. (|26|I - H27|I differ from 
the definition Eq. 123|l by the numerical coefficients in 
the brackets, since the angle integrals are different in 
each regime. However, for rough estimates, one can use 
Eq. 



IV. TWO-PAIR ABSORPTION 

Another mechanism effective for the absorption of plas- 
mons with energies {ujq <C Ep) is associated with cre- 
ation of two electron- hole excitations^. The probability 
of absorption or excitation of two electron-hole pairs by 
an electric field Eq. ^ can be found with the help of 
the Golden rule by treating the field strength and the 
electron-electron interaction Vq perturbatively. 



dw = 2T:\C\'6{ip + Ck - Cp' - ^k' + to) 

(fp'(fk' 



x^(p-Fk-p'-k' + q)- 



{2ny 



(28) 



Here p, k and p', k' are, respectively, the initial and final 
momenta of two electrons. The matrix element C de- 
pends of the initial spin state of the two electrons. For a 
singlet state the transition matrix element is 
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-^0,0 - e<po I 7 z. \ H V ■ V , (,^yj 

. 4p - c;p+q + ^ Sp' ~ 4p' -<A - ^ c;k - 4k+q + ^ c;k' " c;k'-q - ^ J 



and for a triplet state it is 



r ^A. \ ^k-k' ~ ^P'-k , Vk~k' - Vp-k' ^P'-P - ^p-k' Vp'-p - Vp'-k I .„„^ 

-t-i.o - e</)o 7 ; \ y 7 J H ■ h . 

,4p-4p+q+^ ?p'-4p'-q-t^ i;k-4k+q+tt) i;k' - 4k'-q - ' ■ 



The first subscript of the matrix element denotes the 
value of the total spin of two electrons, while the second 
subscript stands for its projection onto the z-direction. 
The matrix element in the triplet channel is the same 
for all three spin states, L\fi = >Ci.±i, since the electron- 
electron interaction is spin-independent. The structure 
of the matrix elements Eqs. H29|l and (|30|l is clear from 
the graphic representation shown in Fig. 3. Again, as in 
Fig. 2, the internal lines correspond to the virtual states 
whose propagators are given by {E— Ey)~^ ^ with E being 
the total energy of the particles and Ey is the energy of 
a virtual state. The total probability of the field absorp- 
tion is given by the sum over the initial and final states 
according to 



1 f d?pcPk 
4 J (,27rj* 

(31) 

Here the coefficient 1 /4 prevents from double-counting of 
the initial and final states. The above formulas are valid 



p+q 



p-q/ p 



k 


k-k , , 
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k-k' ^ 
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P' 
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P' 
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P-P' ^ 
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P-P ^ 



k+q 



k-q 



FIG. 3: Graphic representation for the two-pair absorption 
amplitude £, see Eqs. (I291 - 13U1 . The insertion of the wavy 
line (external wave or plasmon) into the two-electron scatter- 
ing process is possible in the shown hero four different ways, 
corresponding to four fractions in Eqs. H29^ - I|3()j . In addition, 
an interchange p' <-> k' in the final state is possible, bringing 
two terms in each of the numerators in Eqs. I29ll - I3mi . The 
relative signs of these terms depend on the spin state of the 
two electrons. 



for any values of uj and q. We now make use of the fact 



that Lu ^ qvp in the region of interest, and obtain 
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•^0,0 = 
•^1.0 



e<j>o 



{A -A c 



(32) 



Here we introduced the following notations: 



-4 = q • [(Q + q)VQ+q - QVq] 
Ax = q- [(p'-k)V-k + (k'- 



p)^^P-^ 



(33) 



and accounted for the momentum conservation evident 
from Eq. I|28|l . The momentum transferred in a collision 
is denoted by Q = p — p'. The term is obtained from 
A by interchanging the momenta p' <~> k', and originates 
from the indistinguishability of colliding particles. 

The transition rates given by Eq. H28() with the ma- 
trix elements Eq. (|32|) determine the total probability of 
absorption Eq. H31(l . As usual, the total probability of 
emission /g is related to the probability of absorption by 
the detailed balance principle, le = Iae~'^^^ . The en- 
ergy dissipation rate dW/dt can be expressed in terms 
of the probability la according to Eq. (|18|l . The optical 
conductivity is then obtained from its relation to dW/dt, 
represented by Eq. ((T^ . 



a'{uj,q) 



(1 - e-"/^) f d'^pd^kd^p'd'^k' 



rn^oj^q^ 



(27r)5 



(^ - Ax/2) (5(^p + ^k - Cp' - Ck' + t^) 
xupUkil - np>){l - nk')5{^ + k - p' - k' + q). (34) 

The optical conductivity Eq. (|34|) vanishes as 
o' {uj, q) oc q^ in the long-wavelength limit g — > 0. This is 
clear from the small-g properties of the amplitudes A and 
Ax- Indeed, the relation ^ cx at g — > follows directly 
from the first line of Eq. (|33|) , and the relation ^ox oc q^ 
can be obtained from the second line of Eq. H33|) and the 
momentum conservation condition. The absence of ab- 
sorption in the long-wavelength limit, (t'(cij, q = 0) = 0, is 
a manifestation of a more general principle^*': a transla- 
tionally invariant electron liquid does not absorb energy 
from an applied uniform ac field. From a technical stand- 
point, the zero value of a'{uj,q = 0) comes as the result 
of a cancellation of the four matrix elements. Fig. 3, of 
which any element is not small by the external wave vec- 
tor q. Recently, Gornyi and Mirlin^-'^ demonstrated a 
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similar cancellation in the diagrammatic calculation of 
the homogeneous conductivity, while Pustilnik et. al. en- 
countered g^-behavior in the calculation of the far tails 
of the density-density correlation function in the case of 
spinless one-dimensional fermions^^. This cancellation 
was overlooked by Reizer and Vinokur in their calcula- 
tion of the attenuation of two-dimensional plasmone^'^. It 
led to a result large by a factor fc|,/q^, in contradiction 
with the above mentioned general principle. 

The contribution of the exchange term to the 
conductivity Eq. (|34|l can be estimated from Eq. (|33|l . 
This contribution involves large transferred momenta 
|p' — k| ~ kp, whereas the small- momentum domain of 
Q = p — p' is important for the contribution of A. Typi- 
cal values of Q in Eq. (|34|l can be estimated from the con- 
servation of energy, Q ^ max{uj,T)/vp. From Eq. H33(l 
we then obtain Acx/A ^ Vq/Vkp- The amplitude Acx 
can thus be neglected in Eq. H34|l . if the interaction po- 
tential is long-ranged, Vq <^ Vkp ■ This condition is met 
for Coulomb interaction at Ci;,r ^ Ep- In this case the 
scattering amplitudes are the same for singlet and triplet 
states. 

Neglecting the exchange contribution one can conve- 
niently represent the optical conductivity in terms of the 
structure factors for free electrons, 

e^sinhf:^] f cPQ , , 

J sinh[#]sinh[^] ■ 

— C50 



The formula H35() describes energy dissipation due to the 
decay into two electron-hole pairs with energies and mo- 
menta J7 + cj, Q + q and — fi, — Q, respectively. The den- 
sities of the electron-hole pairs are determined by the 
structure factors x" ■ 

In the case of Coulomb interaction, [q- (Q + q) Vq+q — 
q- Q Vq]^ cx q'^/Q^- The momentum integral in Eq. H35() 
is, therefore, singular at small values of Q. If the smallest 
value of Q = uj/2vf, allowed by the energy conservation 
condition, becomes smaller than the inverse screening ra- 
dius K, then the screening of Coulomb interaction must 
be accounted for. This is achieved by the replacement of 
the interaction potential Vq in Eq. (|35|l by its screened 
value. As in the case of the phonon-assisted conductivity, 
it is sufficient here to account for the dielectric function 
by using its static limit, Vq e*VQ/e{0, Q). An empha- 
sis on the small-Q domain comes also from the structure 
of the density of states of the electron-hole pairs. Indeed, 
integration over frequencies Q, in Eq. H35|l results in an 
additional factor c>c . Therefore, if the calculation is 
to be performed with the logarithmic in Kp/u) accuracy, 
we can neglect the difference between Vg+g and Vq, and 



simplify Eq. H35() to 



e2(j2sinh[^] [ QdQ 







£2(0, Q) cosh [P] - cosh [^] ' ^ ' 

where fl± = flzLoj/2. In the long- wavelength limit, Q <C 
kp, the structure factor has the form 



x"{^,Q) 



0{Qvp - \n\). (37) 



To evaluate the conductivity Eq. (|36|l . it is technically 
more convenient to change the order of integrations and 
evaluate first the integral over the momentum. With the 
logarithmic accuracy, the integral is cut at the upper limit 
Q ~ K, leading to the following integral over frequency. 



cr'{uj,q) = 



e^q^ sinh [^] 



dn 



cosh [^y] - cosh [y] 



■In 



KVp 



(38) 



Utilizing the fact that the logarithm is a slow function of 
its argument, one can evaluate the integral in Eq. H38(l 
with the logarithm assigned its value at the characteristic 
frequency A = max(a;, 27rr). The two-pair contribution 
to the optical conductivity takes the form 



a'{uj,q) 



e-q- 



[c.2 + (2.T)2]ln[^]. (39) 



Finally, the plasmon attenuation is found from its rela- 
tion Eq. H14|l to the optical conductivity. 
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K + (2.T)2] In [. 



KVp 



(40) 



where Ag = max(a'g, 27rr). The expression 1)40(1 weakly 
depends on the strength of the electron-electron interac- 
tion. This is due to the fact that the principal contri- 
bution to the conductivity comes, with logarithmic accu- 
racy, from the momentum range Q < n, where interac- 
tion is screened. Our calculation is, therefore, correct as 
long as the Fermi momentum exceeds the inverse screen- 
ing radius, kp ^ k. In a realistic 2DEG, however, the 
two quantities are of the same order of magnitude. In 
this general case it is no longer a good approximation to 
disregard the scattering with large momenta Q ^ kp and 
to neglect the exchange effects represented by the ampli- 
tude Acx in Eq. (|34|l . One should expect the numeri- 
cal coefficient in Eqs. (|39|l - (|40|) to change and become 
interaction-dependent for k}^ kp. 
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V. NUMERICAL ESTIMATES FOR A GaAs 
HETEROSTRUCTURE 

In this Section we compare the effectiveness of different 
mechanisms contributing to the attenuation of plasmons. 
We concentrate on the case of low temperature, and con- 
sider attenuation 7^ as a function of the plasmon energy 
hujq. 

In the important case of a GaAs heterostructure, the 
material parameters are^iSii^: deformation potential cou- 
pling constant S = 2.2 x 10~^^ J, piezoelectric con- 
stant /ii4 = 1.38 X 10^ V m~^, longitudinal sound ve- 

s~^, transverse sound velocity 



locity si ~ 5.2 x 10'^ m 1 
St ~ 3.0 X 10^ m s^^, crystal density p = 5.3 x 10^ 
kg m"'^, and the dielectric constant e* = 12.8. For esti- 
mates, we take a typical Fermi momentum, fcj? « 1 x 10^ 
m^^, which corresponds to the Fermi energy Ep « 5.6 
meV. 

Using Eqs. © and (|23|) we find the dimensionless 
electron-phonon coupling strengths: rjd = 0.03 for the 
deformation potential, and rjp = 0.09 for the piezoelec- 
tric coupling. The characteristic frequencies ujg and cj^i 
introduced after Eqs. ||HJ) and H22|l respectively, are of the 
same order, Tit^K ~ 0.7 meV, and TilDs ~ 0.3 meV. 

Comparison of Eqs. lO, and with each 

other shows that the plasmon-phonon conversion mecha- 
nism yielding 7^ cx w^, see Eq. dominates the attenu- 
ation only at very low plasmon energies, fiLOq ^ 0.02 meV. 
At higher energies, the phonon- assisted Landau damping 
via piezoelectric coupling prevails, and 7^ 



as tkUq Sj 2ti>. 



as long 

0.7 meV. At higher energies, the piezo- 
electric contribution saturates, see Eq. H2()|l . while the 
contribution coming from the deformation potential in- 
teraction, Eq. 124|l . monotonically increases, 7^ o 
It becomes the dominant one at huq ^ 1 meV, 
exceeds the two-pair mechanism, Eq. (|40|) . by a factor 
~ {EF/oJqY ill the entire range hujq < Ep- We are not 
considering higher plasmon energies, where the conven- 
tional Landau damping significantly contributes to the 
plasmon line width. 



and 



VI. CONCLUSIONS 

We presented a theory for the real part of the opti- 
cal conductivity cr'(w, q) of a ballistic 2DEG at finite 
wave-vectors and utilized it for calculation of a plas- 
mon width. Collisionless energy dissipation (Landau 
damping) in a degenerate plasma is exponentially small 
~ exp [— TOaj^/2g^r] in the low-temperature and long- 
wavelength limit. That prompted us to consider a num- 
ber of mechanisms which would be sub-leading, should 
conservation laws allow for the Landau damping. These 
mechanisms are the plasmon-phonon conversion, see Sec- 
tion II, phonon-assisted creation of an electron-hole pair, 
see Section III, and two-pair absorption, see Section IV. 
In the domain of low plasmon energies and low temper- 
atures, hujq,T ^ Ep, all the found contributions have 



power-law (in ojq and T) asymptotes, see Eqs. ((S)), (|21ll . 
(|12J), (Unil, (113, and The smallest exponent be- 

longs to the conversion mechanism, which therefore dom- 
inates the plasmon line width at the lowest temperatures. 
(Wc do not discuss here a competing contribution from 
the impurity scattering.) The relative importance of the 
other considered mechanisms at higher temperatures de- 
pends on the material parameters. In the important case 
of a GaAs heterostructure, apparently the two-pair mech- 
anism does not become the leading one, see Section V. 
This mechanism, however, results in a "universal" con- 
tribution, see Eq. H40|) , which is only weakly (logarithmi- 
cally) dependent on the parameters of the host material. 
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APPENDIX A: PLASMON-PHONON 
CONVERSION RATE IN THE GOLDEN RULE 
APPROXIMATION 

The plasmon-phonon conversion rate, Eq. |(7J), can also 
be derived from the Golden rule formalism similar to the 
phonon-assisted Landau damping and the two-pair atten- 
uation rate. The probability of electric field absorption 
accompanied by the phonon emission. Fig. 4, is given by 
the expression, 

fk) 

dwxk = 27re^|0o|\^(^,g) ^ S{u; - cjxk) 

x%-k)|) — . (Al) 
After summation over phonon states, the total probabil- 




FIG. 4: Transition amplitude for the direct conversion of the 
electric field (wavy line) to a phonon (zigzag line). The loop 
stands for the electron polarization operator xi^tl)- 



ity of the field absorption becomes 



(A2) 



8 



The total probability of the field emission (accompa- 
nied by the phonon annihilation), /g, is obtained from 
Eq. I)A2|I by changing 1 + N\k Nxk- As always, the 
two probabilities are related by the detailed balance con- 
dition, /e — /aC""/"^. The corresponding contribution to 
the optical conductivity is found from Eqs. H18|l and 



Finally, the plasmon attenuation is determined by 
Eq. ini). 



7. = %^E-a'^a(«'V«>)' (A4) 



4p 



cr'{uj,q) 



E 



(A3) which is equivalent to Eq. 0. 
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